Mon. Not. R. Astron. Soc. 000, [THS] (2001) Printed 1 February 2008 (MN WI^ style file v2.2) 



On the width of cold fronts in clusters of galaxies due to 
conduction 



o 
o 



F. Xiang^, E. Churazov^'^, K. Dolag^, V. Springel\ A. Vikhlinin^ 

^ Max- Planck- Institut fiir Astrophysik, Karl-Schwarzschild-Strasse 1, 85741 Garching, Germany 

^ Space Research Institute (IKI), Profsoyuznaya 84/32, Moscow 117810, Russia 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 



1 February 2008 



o 



> 
in 
in 
o 
o 

\6 
o 

o 



ABSTRACT 

We consider the impact of thermal conduction in clusters of galaxies on the (unmag- 
netizcd) interface between a cold gaseous cloud and a hotter gas flowing over the cloud 
(the so-called cold front) . We argue that near the stagnation point of the flow conduc- 
tion creates a spatially extended layer of constant thickness A, where A is of order 
^ ^ kR/U , and R is the curvature radius of the cloud, U is the velocity of the flow at 
infinity, and k is the conductivity of the gas. For typical parameters of the observed 
fronts, one finds A <C i?. The formation time of such a layer is ~ R/U . Once the layer 
is formed, its thickness only slowly varies with time and the quasi-steady layer may 
persist for many characteristic time scales. Based on these simple arguments one can 
use the observed width of the cold fronts in galaxy clusters to constrain the effective 
thermal conductivity of the intra-cluster medium. 
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1 INTRODUCTION 

Chandra observations of galaxy clusters often show sharp 
discontinuities in the surface brightness of the hot intra- 
cluster medium (ICM) emission (Markevitch et al., 2000, 
Vikhlinin, Markevitch, Murray, 2001, see Markevitch & 
Vikhlinin 2007 for a review). Most of these structures have 
lower temperature gas on the brighter (higher density) side 
of the discontinuity, contrary to the expectation for non- 
radiative shocks in the ICM. Within the measurement uncer- 
tainties, the pressure is continuous across these structures, 
suggesting that they are contact discontinuities rather than 
shocks. In the literature these structures are now called "cold 
fronts" . 

There are several plausible mechanisms responsible for 
the formation of such cold fronts, all of them involving rel- 
ative motion of the cold and hot gases. Below we will con- 
sider the case of a hot gas flow over a colder gravitationally 
bound gas cloud, which is a prototypical model of a cold 
front. In such a situation one expects that ram pressure of 
the hotter gas strips the outer layers of the colder cloud, 
exposing denser gas layers and forming a cold front near the 
stagnation point of the hot flow (Markevitch et al., 2000, 
Vikhlinin et al, 20Qla, Bialek, Evrard and Mohr, 2002, Na- 
gai & Kravtsov, 2003, Acreman et al., 2003, Heinz et al., 
2003, Asai, Fukuda & Matsumoto, 2004, 2007, Mathis et al.. 



2005, Tittley & Henriksen, 2005, Takizawa, 2005, Ascasibar 
& Markevitch 2006) . 

Some of the observed cold fronts are remarkably thin. 
For example, the width of the front in Abell 3667 (Vikhlinin 
et al., 2001a) is less than 5kpc, which is comparable to the 
electron mean free path. Given that the temperature changes 
across the front by a factor of ~ 2, thermal conduction (if not 
suppressed) should strongly affect the structure of the front 
(e.g. Ettori & Fabian, 2000). In fact, suppression of conduc- 
tion by magnetic fields is likely to happen along the cold 
front since gas motions on both sides of the interface may 
produce preferentially tangential magnetic field, effectively 
shutting down the heat flux across the front (e.g. Vikhlinin 
et al., 2001b, Narayan & Medvedev, 2001, Asai et al., 2004, 
2005, 2007, Lyutikov 2006). While magnetic fields are hence 
likely to play an important role in shaping cold fronts, it is 
still interesting to consider the expected structure of a cold 
front in the idealized case of an unmagnetized plasma. 

The structure of this paper is as follows. In Section 2, 
basic equations are listed and a toy model of a thermally 
broadened interface between cool and hot gas is discussed. 
In Section 3, we present the results of numerical simulations 
of hot gas flowing past a cooler gas cloud. In Section 4, 
we discuss how limits on the effective conductivity can be 
obtained for the observed cold fronts. Finally, we summarize 
our findings in Section 5. 
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2 THERMAL CONDUCTION NEAR THE 
STAGNATION POINT OF THE FLOW 

2.1 Basic equations 

We parameterize the isotropic tliermal conductivity k as 



k = f X ko, 



(1) 



where / < 1 is the suppression coefficient of the conductivity 
relative to the conductivity ko of an unmagnetized plasma 
(Spitzer 1962, Braginskii 1965): 



ko = 4.6x10' 



(iQSk) 



/lnA\-i 
I 40 J 



erg cm ^ s ' K ', (2) 



where T is the gas temperature, and InA is the Coulomb 
logarithm. 

If the scale length of temperature gradients is much 
larger than the particle mean free path, then saturation of 
the heat flux (Cowie & McKee, 1977) can be neglected and 
the evolution of the temperature distribution can be ob- 
tained by solving the mass, momentum and energy conser- 
vation equations with the heat diffusion term V • fcVT in 
the energy equation (e.g. Landau & Lifshitz, 1959): 
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where p is the gas density, p is the gas pressure, g is the grav- 
itational acceleration, and v is the gas velocity. We adopt an 

ideal gas with 7 = 5/3, where e = 

and p = -E—kBT. 

In the next section we first consider the simplified case 
of passive scalar diffusion in a time independent velocity 
flow, while in Section |3] we discuss numerical solutions of 
the above equations. 
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2.2 Toy model 

Churazov & Inogamov (2004) noted that the behaviour of a 
conducting layer in cold fronts should be similar to the be- 
haviour of a viscous layer near a plate or near the surface of 
a blunt body (see e.g. Batchelor, 1967). When the fluid is ad- 
vected along the surface, the thickness of the layer grows in 
proportion to the square root of the advection time. Near the 
stagnation point, the velocity of the flow increases linearly 
with the distance from the stagnation point and the charac- 
teristic advection time is approximately constant. Therefore 
the thickness of the layer can also be approximately con- 
stant. Below we provide a more rigorous justification of this 
picture. 

Let us consider the simple case of diffusion of a passive 
scalar 7/; in a potential flow of an incompressible fluid. The 
diffusion coefficient D is assumed to be constanlQ and the 



' We use the notation D in this section for constant diffusion 
coefficient to distinguish it from the temperature dependent heat 
conductivity k. 



velocity field is known and constant with time. The diffusion 
equation 



9^ ^ 
— + V 
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{vip) - DAtp = 



(6) 



is supplemented by static boundary conditions at the surface 
of the body and at large distance from the body. For a steady 
state solution (^ = 0) and for an incompressible fluid (V ■ 
V = 0) the above equation reduces to 

w ■ - DA^ = 0. (7) 

In the simplest case of a uniform flow along the "heated" 
plate (Fig.[T]left), Vx = u = const and Vy = 0. At sufficiently 
large distance from the leading edge of the plate, the deriva- 
can be neglected and equation ((T)) can be written 
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An obvious solution in the form ip — f{y/y/x) is given by 



i/i — {ipi — 1P2) Erf 



2Dx- 
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(9) 



where ■02 and i/'i a-re the values of the scalar at the plate 
and at infinity, re spect ively. The width of the interface is 
therefore Ay = \f^^-x and it increases with the distance x 
from the leading edge of the plate as ^/x. Since it takes a 
time t = x/u for the gas to flow from the edge of the plate to 
a given posi tion x, the width of the diffusive layer is simply 
~ ^/m = ^J~DxJu. 

Consider now a potential flow into a 90 degrees cor- 
ner (Fig. [1] middle), governed by the velocity potential 
= Ar^ cos 2Q (see e.g. Lamb 1932, for various examples 
of potential flows). Here r is the distance from the corner 
and Q is angle from the horizontal axis. In this case the veloc- 
ity components are = 2 Ax and Vy — —2 Ay. An obvious 
solution to equation (0 is then 



tp = {ijji — 1^2) Erf 




(10) 



with the width y = of the interface being independent 
of X. The reason for this behaviour is clear: the accelera- 
tion of the (incompressible) fluid along the interface causes 
a contraction of the fluid elements perpendicular to the di- 
rection of the acceleration. While diffusion is trying to make 
the interface broader, the motion of the fluid towards the 
interface compensates for the broadening of the interface, 
and a steady state is reached (Fig. [1] middle). 

The potential flow past a cylinder or sphere behaves 
qualitatively similar (Fig. [T] right). Indeed, in the vicinity of 
the stagnation point (for 6 <^ 1), the radial and tangential 
components can be written as (flow is from the right to the 
left, angle is counted clockwise from the —x direction): 
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for a cylinder and 
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Figure 1. Examples of simple potential flows: uniform flow past a plate (left), flow into a 90 degree corner (middle) and flow past 
a sphere (right). Streamlines are shown as thin solid lines. The shaded areas schematically indicate the width of the layer formed by 
diffusion. 
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for a sphere. Here U is the velocity at infinity, R is the radius 
of the cylinder or sphere, rj — r — R and — RsinO. 

In the same approximation as for the cases discussed 
above (where the spatial derivative of -0 along ( is neglected) 
the diffusion equation reduces to 



dip 
or 



(13) 



and the width of the interface over the radius is set by the 
diffusion coefficient D and the coefficient C in the relation 
Vr = —Ct], yielding 



cylinder 



(14) 



y^D^ sphere. 



In this case the width of the interface is also constant along 
the surface of the cylinder or sphere (Fig. [T] right). 

One can also consider a closer analogue of a flow past a 
spherical cloud by extending the solution for a potential flow 
into the inner part of the cylinder or sphere, as illustrated in 
Fig. [T] In this model there is a circulation flow of gas inside 
the cloud, and the tangential component of the velocity is 
continuous across the boundary while the normal component 
is zero at the boundary. We can further allow for different 
densities pi and p2 outside and inside of the boundary if 
aU velocities inside are scaled by a factor p\l p2- The re- 
sulting configuration can be considered as an idealized (and 
unstable) analogue of a hot flow past a colder cloud in the 
absence of gravity (see also Heinz et al., 2003). Allowing dif- 
ferent diffusion coefficients D\ and D2 in the ffow outside 
and inside the boundary, and requiring the solution ■(/) and 
its spatial derivative to be continuous across the interface, 
yields the following solution in the vicinity of the stagnation 
point: 



ipm) Erf ( \\ ^^{t - '"0) 1 + V*™, outside. 
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Here ro is the radius of the boundary , an d "02 are the val- 
ues far from the interface, C2 — Ciy^ P\l P2, and Ci = 2U /R 
for a cylinder or Ci = 3U/R for a sphere, respectively. The 
width of the interface is again constant along the boundary. 

The same answer is obviously valid for any idealized 
flow of this type: near the stagnation point the width of 
the "heated" layer does not change along the surface of the 
body. Real cold fronts are of course much more complicated 
structures. However, the acceleration of the flow along the 
interface and the simultaneous contraction in the perpen- 
dicular direction are generically present also here. It can 
therefore be expected that the width of the interface will 
be similarly constant in real cold fronts. A simple extension 
of the above toy model can be obtained by allowing for gas 
compressibility and a temperature dependent conductivity, 
i.e. by considering the full system of equations ©-((SJ with 
conductivity according to eq. ([2]). An expansion of heated 
layers and simultaneous contraction of cooled layers on the 
other side of the interface will certainly modify the flow, but 
for the transonic flows of interest here we might expect that 
the results obtained for a toy model will still be approxi- 
mately valid. In the next section we verify this prediction 
using numerical simulations. 



3 NUMERICAL SIMULATIONS 

For our numerical experiments, we used the TreeSPH code 
GADGET-2 (Springel, 2005) combined with the implementa- 
tion of thermal conduction discussed by Jubelgas, Springel 
& Dolag (2004), which accounts both for the saturated and 
unsaturated regimes of the heat flux. 

The simulations were intended to illustrate a simple toy 
model, described in section \%2\ rather than to provide a re- 
alistic description of the observed cold fronts. The specific 
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goal was to see the impact of the flow stretching near the 
stagnation point on the width of the interface set by conduc- 
tion. With this in mind we intentionally restricted ourselves 
to a 2D geometry and an unmagnetized plasma. For a 3D 
calculation of magnetized clouds see Asai et al. (2007). The 
self gravity of gas particles was also neglected in our ide- 
alized simulations and all gas motions were happening in 
a static gravitational potential. Given that the typical gas 
mass fraction in clusters is of order 10-15 per cent, the self 
gravity of gas particles is likely to be a second order ef- 
fect. A more significant simplification is the assumption of 
a static potential, since at least some of the cold fronts are 
caused by cluster mergers where strong changes of the poten- 
tial are possible. Formation of cold fronts in the appropriate 
cosmological conditions was considered by e.g. Bialek et al. 
(2002), Nagai & Kravtsov (2003), Mathis et al. (2005), see 
also Tittley & Henriksen (2005) and Ascasibar & Markevitch 
(2006). Our illustrative 2D simulations, described below, can 
be viewed as a "minimal" configuration which allows us to 
see the effect of flow stretching and to extend the toy model 
to the case of a compressible gas and a temperature depen- 
dent diffusion coefficient. 



of initial particle positions and it does not mean that the 
number of particles is insufficient to properly resolve the 
cold front. A comparison of three runs with 3 10^, 10^ and 
3 10*" particles with the same initial condictions and the con- 
duction suppression coefficient / — 0.01 yielded practically 
undistinguishable results in terms of the cold front struc- 
ture. As these ripples do not affect the overall structure of 
the flow no attempt was made to correct the initial condi- 
tions for this effect. These ripples are also a useful visual 
indicator of the impact of thermal conduction on the small- 
scale temperature structures in the flow (see Fig. [Sjl . 

Our choice of initial conditions has been motivated by 
the cold front in the cluster Abell 3667 (Vikhlinin et al., 
2001), but we did not try to accurately reproduce all the 
observed properties of this cluster. In particular, the location 
and the strength of the shock in our model need not be the 
same as in A3667. Nevertheless, the most important feature 
of A3667 - a cool gas cloud inside a hotter and less dense 
flow - is present in our simulations, allowing us to study 
the impact of thermal conduction on the interface between 
the cloud and the flow. To this end we carried out four runs 
where the coefficient of the thermal conduction efficiency 
was set to / = 0, 0.01, 0.1, and 0.5, respectively. 



3.1 Initial conditions 

Our 2D simulations of cold fronts in clusters were carried 
out in a 8x4 Mpc periodic box. We represented the cluster 
with a static King gravitational potential of the form 



In [a; -f Vl + a;^ 
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with a = 810 km s 



(15) 

■ r/vc and Tc = 300 kpc. The initial 



temperature and density distributions were set to 
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and kTi = nrripa w 4 keV. Thus the temperature and 
density make a jump at Xout , while the pressure is contin- 
uous. In our runs, T2 — S keV, a;out = 1, /i = 0.61, and 
pi — 6.6 X 10^^^ g cm^"^. The gas velocity was set to zero 
for X < Xout and to u = 2000 km s~^ for x > Xout. The 
corresponding Mach number relative to the hot 8 keV gas 
is ~1.3 (neglecting further acceleration of the flow in the 
cluster potential). 

10^ gas particles of equal mass were distributed over the 
computational volume as a Poissonian sample of the initial 
density distribution described above. The gas temperature 
for each particle was set to Ti or T2 depending on the posi- 
tion of the particle. The Poissonian noise introduced by this 
procedure leads to small-scale (and small amplitude) pres- 
sure/entropy perturbations in the initial conditions. Shortly 
after the beginning of the simulations, the over-dense regions 
expand and create a pattern of ripples in the temperature 
distribution (see Fig. [21 top panel). We stress that the pres- 
ence of these ripples is a direct consequence of the choice 



3.2 Results 

In Fig. [2] we show the temperature distributions for all four 
runs 1.3 Gyr after the start of the simulations. Immediately 
after the beginning of the simulations a shock starts to prop- 
agate upstream through the hot flow, forming a clearly visi- 
ble bow shock. Because of the acceleration in the cluster po- 
tential, the Mach number of the shock is ~1.7 (rather than 
1.3) and the temperature behind the shock is also rather 
high (10-15 keV). 

The cool cloud is flrst pushed back by the ram pressure 
of the gas and then (slowly) oscillates near an equilibrium 
position. At 1.3 Gyr, there are still some residual motions 
clearly associated with the specific initial conditions, but 
these motions are quite gradual. This can also be seen in 
the gas velocity field, which is plotted in Fig. [3] It shows a 
clearly visible velocity jump at the shock front, and inside 
and around the cloud, circular motions are present, broadly 
resembling the velocity field shown in Fig. [T] Such circular 
motions inside the cloud lead to a transport of low entropy 
gas from the centre of the cloud towards the stagnation point 
(Heinz et al., 2003). As a result of adiabatic expansion of the 
transported gas its temperature drops to ~ 3 keV, below the 
initial value of 4 keV. 

Since after 1.3 Gyrs much of the relaxation from the ini- 
tial state already took place, we compare the runs with dif- 
ferent conductivity at this time. The effect of increasing the 
efficiency of thermal conduction is clearly visible in the snap- 
shots shown in Fig. (2] First of all, small scale temperature 
variations present in the initial conditions are smoothed out 
in all runs where thermal conduction is present. Secondly, 
with the increase of / the interface separating the cloud and 
the hot flow becomes less and less sharp. This is seen more 
explicitly in the temperature proflles across the interface 
(along the symmetry axis of the cloud), which are shown 
in Fig. m In this flgure (and in the subsequent flgures), the 
distance (plotted along the abscissa axis) is measured from 
the approximate centre of the cloud. Since the cloud is not 
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Figure 3. Snapshot of the temperature distribution at t = 1.3 Gyr after the start of a simulation with conductivity coefficient / = 0.1, 
with the velocity field superposed. 



perfectly spherical, its position varies slowly with time. As 
its centre is hence not accurately known, all profiles shown 
in Fig. |4] were shifted along the abscissa axis to have the 
temperature value 6 keV gas at the same position. 

The sharpest profile corresponds of course to the run 
without conduction, and in this case some small-scale fiuc- 
tuations left over from the initial conditions can still be seen 
in the profile. This run also sets a useful benchmark for com- 
parison with the other simulations, for example, it indicates 
the numerical resolution available for representing the in- 
terface. We see that for values of / larger than 0.01 the 
impact of thermal conduction on the width of the interface 
can be well resolved with our numerical setup. For runs with 
/ — 0.01, 0.1 and 0.5, the small-scale fiuctuations in the tem- 
perature distribution are absent and the effective thickness 
of the interface gradually increases. 

We can now verify our simple predictions based on the 
toy model of diffusion of a passive scalar in a potential flow. 
We first consider our finding that after an initial settling 
time of order R/U the interface evolves to a quasi-steady 
state. This is indeed seen in Fig. (6] where the tempera- 
ture profiles along the symmetry axis are shown for t = 0.2, 
0.4, 0.8, and 1.3 Gyr since the beginning of the simulations. 
While there is clear evolution of the profile (e.g. in terms 
of the maximal or minimal temperatures) the shape of the 
interface is very similar at all times. 

Another expectation is that the thickness of the inter- 
face is the same along the interface (as long as the distance 
from the stagnation point is much less the the cloud curva- 
ture radius). Indeed, the profiles measured at different dis- 
tances from the stagnation point (Fig. [T]) look very similar. 
In this figure, the profiles were calculated along directions 
making different angle with respect to the symmetry axis of 
the cloud (red lines in Fig. [2]). This means that when deriv- 
ing an effective conduction coefficient from the observed cold 
fronts one can use the profile averaged over a large part of 
the interface, rather than being constrained to small sectors 
of the front. 

Thus the results of the numerical simulations are 
broadly consistent with the expectations derived earlier: 



once the front is formed, it has a width constant in time 
and constant along the interface. 



4 SIMPLE ESTIMATES OF THE INTERFACE 
WIDTH 

We are now looking for a simple ID problem for a compress- 
ible fiuid which has a solution that can provide a qualitative 
approximation of the cold front structure. The analogy with 
the problems considered in Section 12.21 suggests that the 
width of the interface should scale as \/Dts, where D is the 
effective diffusion coefficient (thermal conductivity) and ts 
is the effective time scale. For the flows in Section 12.21 a 
reasonable choice was ts ~ R/U, where R is the curvature 
radius of the interface and U is the velocity of the flow at 
infinity. Of course, this result was derived for a potential 
flow, and for more realistic cases it might be more correct 
to recast ts in the form ts — (^) , where Vr is the veloc- 
ity component perpendicular to the interface. Indeed, from 
Eqs. (|13|) and (|14|l it is clear that that this quantity (i.e. the 
gradient of the radial velocity component) enters the expres- 
sion of the interface width. We can hence try to obtain an 
approximate solution for the front structure by considering 
a ID time dependent diffusion equatior0, starting from a 
Heaviside step function for the initial temperature distribu- 
tion and taking the solution at time ts- The hope is that for 
this choice of ts the most basic properties of the interface 
structure will be captured. 

We also assume that all velocities in the vicinity of 
the interface are small compared to the sound speed, all 
quadratic terms in v can be neglected and that the pressure 
is approximately constant across the interface. Then the gas 
density is p = po , where poTo = Po is fixed by the initial 
pressure. In this approximation, the heat diffusion equation 



^ An alternative approach is to incorporate the velocity field of 
the base flow into the energy conservation equation and to look 
for a steady state solution. 
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Figure 4. Temperature profiles along the symmetry axis of the 
cloud at t = 1.3 Gyr after the start of the simulations, for different 
conduction suppression coefficients equal to / = 0, 0.01, 0.1, and 
0.5. The distance is measured from the centre of the cold cloud. 




Figure 2. Snapshots of the temperature distribution for a 2D 
flow of hot (8 keV) gas past a cold (4 keV) cloud. The image 
sizes are 6.7 by 2.7 Mpc. The top panel corresponds to a moment 
shortly after the start of the simulations. In the other panels, the 
temperature distributions at t = 1.3 Gyr after the beginning of 
the simulations are shown as a function of the strength of ther- 
mal conduction. The conduction suppression coefficient for these 
panels is / = 0, 0.01, 0.1, and 0.5, respectively. Thus the second 
panel shows the run without conduction, while the bottom panel 
corresponds to a conductivity equal to half the Spitzcr-Braginskii 
value. The temperature structure of the interface plotted in sub- 
sequent figures was measured along the red lines shown in one of 
the panels. 



reduces to 



dT 



dt dx dx 



4 m 



(18) 



where A — - — ^-r^. This equation is very similar to the 



standard diffusion equation in solids, except for the second 
term in the r.h.s. which accounts for gradual expansion of 
the heated gas, and for the contraction of the cooled gas in 
order to maintain constant pressure across the interface. 
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Figure 5. Gas density profiles along the symmetry axis of the 
cloud at t = 1.3 Gyr after the start of the simulations, for different 
conduction suppression coefficients equal to / = 0, 0.01, 0.1, and 
0.5. The distance is measured from the centre of the cold cloud. 
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Figure 6. Temperature profiles along the symmetry axis of the 
cloud after t = 0.2, 0.4, 0.8 and 1.3 Gyr since the start of the simu- 
lations. The conduction suppression coefficient was set to / = 0.1. 
The distance is measured from the centre of the cloud. In order to 
compensate for the gradual changes in the interface shape, each 
profile was shifted along the X-axis so that the rising part of the 
temperature profile has the same abscissa for all curves. 
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Figure 7. Temperature profiles at 1.3 Gyr for / = 0.1 along 
directions making an angle of 0, 15, 30, or 45 degrees with respect 
to the symmetry axis of the cloud (see also Fig.[2]|. 



Equation (|18|l can be readily integrated. We set the ini- 
tial values of the temperature to 3 and 15 keV on the two 
sides of the interface, respectively, and the electron density 
to 10~^ cm"'' on the cool side to approximately reproduce 
the properties of the simulated front, as shown in Fig. [7] 

In Figure El we plot the solution of the equation at times 
t = 0.1, 0.4, 0.9, and 2.0 Gyr for / = 0.1 together with the 
temperature profile obtained in the SPH simulations. The 
results of the numerical simulations best correspond to the 
solution of Eqn. (|18|) for t ~ 0.5 — 0.7 Gyr. For comparison, 
the ratio of the cloud radius to the flow velocity at infinity is 
R/U ~ 1.5 X 10* yr, while ts evaluated from the velocity pro- 
file obtained in the simulations is ts = (^j^) ~ 5 x 10* yr. 
The difference in the estimated width of the front based on 
the R/U ratio compared with a more detailed treatment of 
the velocity field is of order 2. This discrepancy (for our 
numerical setup) is largely caused by i) the drop of the ve- 
locity at the shock and ii) differences between the veloc- 
ity field obtained in the simulations and that in potential 
flows of incompressible fluids, as considered in Section 12.21 
Of course, Ean. (|18| l by itself is only a crude approximation of 
the problem. Nevertheless, even our simplest estimate pre- 
dicts the width of the interface within a factor of 2 of the 
value derived from direct numerical simulations. 

Note also that because of the expansion of the heated 
gas, the actual contact discontinuity does not necessarily ex- 
actly coincide with the "visible" boundary of the cool cloud. 
Indeed, colder gas of the cloud expands while the hotter am- 
bient gas contracts, shifting the discontinuity away from the 
cloud centre. At the same time, the sharpest edge will be 
observed in places where the temperature gradient is large. 
If the gas on the two sides of the contact discontinuity has 
different abundances of heavy elements then the true posi- 
tion of the contact discontinuity can be determined from the 
abundance gradient. This exercise however requires data of 
very high quality. 

We thus flnd that an order of magnitude estimate of the 
width of a spherical cold front can be written as 
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(19) 



where the factor S accounts for all departures introduced 
by the approximations involved in our simplest model con- 
sidered in Section 12.21 Based on our numerical simulations 
we have 5 ~ 0.5. Plugging in fiducial values for the A3667 
cluster (Vikhlinin et al., 2001) one gets: 

5 



Ar 40 
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(20) 



and 



2x 10"3 cm^3 ■ 

This value is a factor of ~ 8 larger than the upper limit 
for the interface width derived from Chandra observations. If 
this discrepancy is solely caused by a suppression of thermal 
conduction, then the factor / has to be less than 1.5 x 10~ . 
Note that the mean free path is A ~ 4 Ts A'^qqj kpc (e.g. 
Sarazin, 1986), which is smaller than the interface width 
evaluated for / = 1. If the effective mean free path of elec- 
trons scales linearly with f = k/ko (while the interface width 
Ar oc ^/f) then for all / < 1 the width of the interface will 
remain larger than the mean free path. Therefore our as- 
sumption of unsaturated heat flux remains valid. For the 
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Figure 8. Comparison of the temperature profile along the sym- 
metry axis of the cloud (for / = 0.1) obtained in numerical simu- 
lations (histogram) and obtained from equation (I18t for ts = 0.1, 
0.4, 0.9, and 2.0 Gyr. 

observed cold front in A3667 the effective mean free path of 
electrons should therefore be factor / = 1.5 x 10~^ smaller 
than in unmagnetized plasma. 

However, we caution that magnetic fields are likely play- 
ing a role in the structure of the interface as suggested by 
theoretical arguments (e.g. Vikhlinin et al., 2001, Lyutikov 
2006) and numerical simulations (e.g. Asai 2004, 2005). Par- 
ticularly important is the stretching of the field lines along 
the interface, which can suppress heat conduction across the 
front or even affect the hydrodynamical stability of the in- 
terface. Note that the heat conductivity depends strongly on 
the topology of the magnetic field since the electron Larmor 
radius is some 10 orders of magnitude smaller than the char- 
acteristic length scales of the problem for typical magnetic 
field strengths at the micro-Gauss level. 

We note that the interface developing in our SPH simu- 
lations may also be stabilized to some extent against small- 
scale fluid instabilities by numerical effects. Across strong 
density discontinuities, SPH has been found to produce spu- 
rious pressure forces that may suppress small wavelength 
Kelvin-Helmholtz instabilities (Agertz et al. 2006). This ef- 
fect is equivalent to a small surface tension and mimics the 
stabilizing influence expected from an ordered magnetic field 
across the front. However, better numerical resolution weak- 
ens this effect and should allow ever smaller wavelengths to 
grow. 

5 CONCLUSIONS 

We have shown that in the presence of thermal conduction 
the width of the interface separating hot gas flowing past a 



cooler gas cloud (a "cold front" in clusters of galaxies) can 
be estimated from the size R of the cloud, the velocity U of 
the gas and the effective thermal conductivity. The structure 
of the interface is established over a period of time ^ R/U, 
while the subsequent evolution is much slower. Moreover, 
the width of the interface is approximately constant along 
the front. We made an illustrative 2D simulations of an un- 
magnetized plasma flow past a colder cloud with gas den- 
sities and temperatures characteristic for the observed cold 
fronts. While being very idealized, the simualtions do show 
that the width remains approximately constant when the 
gas compressiblity and the temperature dependence of the 
conductivity is accounted for. 

This implies that one can use much of the visible part 
of the interface in order to assess the effective thermal con- 
ductivity of the gas. For the cold front in Abell 3667, the es- 
timated width of the interface is 40/'^"'' kpc, where / is the 
conduction suppression coefficient (relative to the Spitzer- 
Braginskii value). This factor / has to be smaller than 0.015 
in order to reproduce the observed limits on the width of 
the interface. This result is consistent with previous sug- 
gestions that magnetic fields play an important role in pro- 
viding thermal isolation of the gases separated by the cold 
front. The idealised description of the interface presented 
here provides a useful method for estimating the effective 
gas conductivity from observations of clusters of galaxies. 
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